Nonadiabatic geometric rotation of electron spin in a quantum dot by 2Pi hyperbolic secant pulses 
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In this paper, the geometric and dynamic phase components of overah phase induced by 2n hyperbolic secant 
pulses in a quantum dot is analyzed. The dependence of two phase components on the ratio of the Rabi frequency 
to the detuning is investigated. Numerical results indicate that only for one resonant pulse the induced overall 
phase is purely the geometric phase. With other values of the ratio the overall phase consists of a nonzero 
dynamic part. The effect of spin precession to decrease the dynamic phase is characterized and discussed by 
analytical and numerical techniques. Utilizing the symmetry relations of the phases, a scheme to eliminate the 
dynamic phase by multipulse control is proposed. By choosing the proper parameter for each pulse, the dynamic 
phases induced by different pulses cancel out. The total pure geometric phase varies from -n to n, which realizes 
the arbitrary geometric rotation of spin. Average fidelity is calculated and the effects of magnetic field and decay 
of the trion state are compared and discussed. The results show the crucial role of weak magnetic field for high 
fidelity (above 99.3%). 

PACS numbers: 78.67.Hc, 03.67.Lx 



I. INTRODUCTION 

Optical control of a single electron spin in a quantum dot 
(QD) is a key ingredient for implementing quantum informa- 
tion processing (QIP) in a scalable solid-state system In 
recent years, great efforts have been made to the physical im- 
plementation of such optical approach. It has been experimen- 
tally demonstrated using oscillating magnetic field generated 
by radio-frequency (RF) pulses on timescales of nanoseconds 
l|2l,i3J. However, ultrafast optical technique is more attractive 
due to its ability of enabling spin rotations to be completed on 
a picosecond timescale which is much shorter than spin coher- 
ent times |4], leading to abundant theoretical proposals [5-13] 
and several experimental achievements il4l4l8ll . 

The proposal T proposed by Economou et at. ifTol [Till is 
based on the attractive feature of the fast laser pulse with a hy- 
perbolic secant temporal envelope 1 19]. After a In sec pulse, 
the two-level system of a QD completes a Rabi oscillation 
and the population returns to its initiate state while having 
acquired an overall phase. This initiate state with the over- 
all phase and a third state which does not couple to the laser 
pulse span the qubit space. The fast spin rotations about the 
optical axis is thus achieved, and the rotation angle is deter- 
mined by the ratio of the Rabi frequency to the detuning, Q/A 
ifToll . In general, the Itt Rabi oscillation may be considered 
as a nonadiabatic cyclic evolution, therefore the overall phase 
consists of a dynamic component and a geometric component 
defined by Aharonov and Anandan |20]. Particularly, when 
the rotation angle (fi - n (f2/A — > oo), i.e. the resonant case, 
the dynamic phase component is zero and the overall phase is 
exactly equals to the geometric phase component. The geo- 
metric phase depends only on some global geometric features 
e.g. the curve in the parameter space, not on the duration of 
interaction. Besides, geometric phases may be robust against 
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dephasing. Although this property is not predicted in general 
systems 02111 . it manifests in particular open systems where 
precession is in the equatorial plane of the Bloch sphere jH^j- 
Therefore the quantum computation based on pure geometric 
phase, so called geometric quantum computing (GQC) ll23ll . 
may have inherent error-resilient advantage. Proposal T is 
thus one of the methods to physically implement GQC. How- 
ever, the cycle-averaged expectation value of the Hamiltonian 
is not always zero with various values of the ratio Q/A, lead- 
ing to a nonzero dynamic component and diverse geometric 
component. In this paper, we will analyze the geometric and 
dynamic phase components induced by the 27r sec pulse in 
detail, through investigating the dependence of the two com- 
ponents on the spin rotation angle and on the ratio Q/A. Two 
methods are utilized: one is to calculate the phases compo- 
nents by substituting the state vector with analytic solutions of 
Ref. 1 10] in intermediate steps, the other is full numerical sim- 
ulation. Besides, the discrepancy of two results is discussed 
and the effect of the other interaction mechanism within the 
system i.e. the spin precession is demonstrated. 

To generalize the proposal T to universal nonadiabatic ge- 
ometric quantum gates, e.g. the arbitrary angle rotation about 
an axis is purely the geometric phase, one needs to avoid or 
remove the dynamic component when (p + tt (Q/A ^ oo). 
In previous studies |24-28], one method is to choose the dark 
states as the cyclic states, thus the dynamic phase component 
is always zero. This scheme has been proposed for nonadia- 
batic GQC with NMR Josephsoii-[unction IH], trapped 
ions |26], and quantum dot systems [27]. The other method 
is the multiloop scheme generalized from the adiabatic evolu- 
tion case ll29i I30II : let the system undergoes evolution along 
several closed loops, thus the dynamic phases accumulated in 
different loops may be canceled, while the geometric phase 
being added. This scheme has been demonstrated in Joseph- 
son junctions and NMR systems [28[. In this paper, we pro- 
pose a scheme to eliminate the dynamic phases in quantum dot 
system. Though being similar with the multiloop method, we 
do not employ the RF pulses to generate the oscillating mag- 
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netic field in our scheme, but employ several picosecond sec 
pulses proposed in the proposal T . The dynamic phases in- 
duced by different pulses cancel out, leading to arbitrary spin 
rotation with pure geometric phase. 

The paper is organized as follows: in Sec. II, we give a brief 
review of the proposal for fast optical rotations of an electron 
spin trapped in a quantum dot, introduce the analytic solutions 
and calculate the expression of the phase components. In Sec. 
Ill, we present numerical results, analyze the phase compo- 
nents and demonstrate the effect of the spin precession to the 
dynamic phase. In Sec. IV, the scheme to eliminate dynamic 
phases by multipulse control is proposed. In Sec. V, we cal- 
culate the fidelity and discuss the effects of magnetic field and 
decay of the trion state. A summary and some prospects are 
provided in Sec. VI. 

II. REVIEW OF THE PROPOSAL AND CALCULATION 
OF THE PHASE COMPONENTS 



0, C7(-oo) = 1, C=(-oo) — 0, the analytic expression of state 




FIG. 1. (Color online) (a) Level diagram of the A system in QD. 
In the Voigt geometry, two lower states are \z) and |z) with electron 
spins in the z direction, coupled by the magnetic field B. The inter- 
mediate state |t) is the trion state with the heavy spin pinned along z 
axis. Circularly polarized light propagates parallel to the QD growth 
direction. Due to selective rules, cr^ pulse only couples |z) and |t>. 
(b) Bloch sphere sketches electron spin and precession of the spin 
around the magnetic field. The clockwise arrows mark the spin pre- 
cession direction. 



The nanostructure employed in the proposal T consists of 
arrays of self-assembled (In,Ga)As/GaAs QDs, each contain- 
ing on average a single electron jsill . In the Voigt geometry, 
the external magnetic field (along the x axis) B - 0.29T Ll7i1 . 
applied perpendicularly to the QD growth direction (z axis) 
which is parallel to the optical axis, as shown in Fig. [Tfb). 
The two lower states \z) and |z) with spins in the z direction 
are coupled through the magnetic field, and they are superpo- 
sitions of spin energy eigenstates in the x direction. When the 
QD absorbs a photon, it's excitated to a trion state which con- 
sists of a singlet pair of electrons and a heavy hole. The hole 
spin is pinned along the growth direction due to strong con- 
finement and spin-orbit interaction. Selection rules determine 
that specific circularly polarized light only couples one spin 
state to one trion state, e.g. cr+ light couples |z) to |r), leading 
to the A-type system shown in Fig.[Tfa). The Hamiltonian of 
system in the interaction picture reads 
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with the basis {\z}, \z}, |t)). lob is spin Larmor precession fre- 
quency and A is the detuning. D.{t) - Q.sech{rit) is the time 
dependant Rabi frequency with a hyperbolic secant temporal 
envelope, where 77 is the pulse bandwidth. 

To obtain analytic solutions, the slow precession approxi- 
mation 1 10] is required: when the Zeeman splitting is much 
smaller than the pulse bandwidth, i.e. ljb rj, the spin 
precession can be neglected within the pulse action. In this 
case the A-type system can be considered a direct sum of 
1 -dimension and 2-dimension systems, with the basis {|z)) 
and Wz), |t)) respectively. By the approach of Rosen and 
Zener [19], the second order equation of probability ampli- 
tudes is transformed to the hypergeometric differential equa- 
tion through change of variable z - (l/2)[tanh(77f) -H 1]. Con- 
sidering the spins are initiahzed in the z direction, Q(-oo) - 



vector is obtained as 



\i^(f)) = F{a,-a;c;z)\z) - -z'' F{a + c,-a + c;\+c;z)\T), (2) 
c 



where F denotes the Gauss hypergeometric series 2F1 1132 
and the parameters a - Q./ri, c = (1/2)(1 -1- iA/rj). 

If and only if Q = 77, the pulse area 2 J Q.sech(rjt) equals 
to In. In this case during the pulse action the system under- 
goes a cyclic evolution between \z} and |r). After the pulse the 
population returns to \z) with an overall phase <p: 



m^)} = e''^\z}, 
Q. 

56 = 2arctan(— ). 



(3) 
(4) 



Considering the qubit space spanned by {\z}, \z}}, the unitary 
operator induced by the 2n sec pulse ITQQ is 



U spin — 



,-/</>/2 Q 

e'«*/2 



(5) 



thus the spin rotation about z axis is achieved. 

Now we demonstrate the calculation of the phase compo- 
nents. The dynamic phase a is the cycle-averaged expecta- 
tion value of the Hamiltonian and the geometric phase y is the 
rest part of removing the dynamic component from the overall 
phase S: 



a=-\ { dtmwm)), 

n J-00 
y— (f> — a. 



(6) 
(7) 



One method (method I) to calculate the phase components is 
to utilize the analytic expression of |i/'(f)) under condition Q = 
rj. We insert Eq. (|2]) into Eq. ^ and rewrite the dynamic phase 
as 
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a = sech2(Qf)[e"''^'(l - tanh(Qr))^(l + tanh(fif))^(A - /Qtanh(Qf)) + c.c.]. (8) 

A2 + J_„ 



In this case the overall phase is calculated following Eq. (|4|. 
We note that the expression of Eq. ([8]) is independent of 
(Ob, resulting from the slow precession approximation. The 
other method (method II) is full numerical simulation, where 
the state vector is calculated through numerical solution of 
Schrodinger equation. For this case the spin precession is 
taken into account. The specific considered parameters are 
demonstrated in Section V. 

The integral interval is infinity in Eq. ©([HI), which is in 
theory for including the whole sec pulse. In experiments the 
stimulating effect of pulses concentrates within the pulse du- 
ration, which is as short as 1.5ps |17[], and the pulse tail on 
a long timescale can be safely neglected. On the other hand, 
especially for the full numerical simulation method where no 
slow precession approximation is made: if the time interval 
is too large, the population transferred to ||) due to spin pre- 
cession is prominent, thus the evolution of system can not be 
considered cyclic and the proposal T would not be ensured. 
Therefore in numerical simulation, we focus on the evolution 
of system during the pulse action. For both methods the inte- 
gral interval is finite, with the same order of magnitude of the 
pulse duration. 



III. ANALYSIS OF THE PHASE COMPONENTS AND 
EFFECT OF THE SPIN PRECESSION TO THE DYNAMIC 
PHASE 

First we analyze the dependence of geometric and dynamic 
phase components on the spin rotation angle (also the overall 
phase) 4>- By the method utilizing the analytic solution intro- 
duced in Eq. (|2]l, numerical results are obtained and showed 
in Fig. |2l The dynamic phase a is a concave function of 0, 
which is symmetrical to the axis of (p = n/l. a is nonzero 
and positive except when (p = 0, n. The maximum of a 
is 2 (blue dot in Fig. |2|i, which is in reasonable agreement 
with the theoretical expression of Eq. ^ equal to 2 when 
<p - 7t/2, Q. - A. The physical origin of the maximum 
dynamic phase is the instantaneous expectation value of the 
Hamiltonian {il/{t)\H\iJ/{t)) attains maximum when system is 
in the eigenstate \ip) - 1/V2(|z) + |t)), which can obtained 
when Q = A. Removing the dynamic component from the 
overall phase leads to the geometric component y, which is a 
convex function of (p. On some interval a > (p, resulting in 
negative values of y. Fig. |2] shows the geometric component 

has a minimum 0.68 when (p k 1.05 (green dot). When 

4> ~ \ .90 the overall phase is equal to the dynamic phase, lead- 
ing to the zero value of geometric component (red dot). The 
concavity of geometric component indicate that only when the 
spin rotated by (p - n, the overall phase is purely the geomet- 
ric phase. For other rotation angles, (p consists of a nonzero 
dynamic component. 

As the review in Sec. II, the rotation angle is determined 
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FIG. 2. (Color online) The dependence of dynamic and geometric 
phase components on the rotation angle (also the overall phase) 0. 
The dashed (blue), dotted (purple) and solid (brown) curves corre- 
spond to the dynamic component a, overall phase ^ and the geomet- 
ric component y, respectively. The blue, red and green dots mark the 
maximum of dynamic phase, zero point and minimum of geometric 
phase respectively. 



by n/A, therefore the phase components are essentially de- 
pendent on Q/A. Fig. [3] shows the log-linear plot of overall 
phase and two phase components as a function of f2/A. The 
dynamic phase component still has a symmetry which is char- 
acterized more clearly. An interesting feature is that O/A and 
its reciprocal leads to the same value of dynamic component, 
e.g. Q'(IO) - Q'(O.l), resulting from the feature of integrand in 
Eq. (|8j. The minimum and zero value of geometric compo- 
nent are obtained when f2/A ^ 0.58 and Q/A » 1.39, respec- 
tively. 
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FIG. 3. (Color online) Log-linear plot of dynamic and geometric 
phase components as a function of the ratio of the Rabi frequency 
to the detuning, H/A. The dashed (blue), dotted (purple) and solid 
(brown) curves correspond to the dynamic component a, overall 
phase (j) and the geometric component y, respectively. 
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Fig.|4]shows the discrepancy of dynamic phase component 
by two methods. The soHd curve corresponds to the result ob- 
tained by method II and the dashed curve by method I which 
is as the same as shown in Fig. |3] Both results are similar 
and in agreement with each other on distribution and variation 
trend. Especially when Q/A <«c 1 two curves basically coin- 
cide. Around Q/A = 1 the discrepancy appears. The dynamic 
phase by method II reaches its maximum earlier than the result 
by method I. With the increase of f2/A, the dynamic phase by 
method II drops faster than the result by method I. But when 
Q/A » 1 two curves approach to each other again. The phys- 
ical origin of the discrepancy is obviously the spin precession, 
which is neglected by method I but considered by method II. 
Due to spin precession about x axis during the pulse action, 
a little population is transferred from \z) to |z). This leads to 
not only degrading the unitarity of the spin rotation operation 
ifTTll . but also the change of system energy distribution. The 
trion state |t) has a higher energy than two lower states |z) and 
\z). The population transferred to |z) will not be transferred 
back to \z) during the pulse action because of ub <^ Q, and ||) 
will not be excitated to the trion state due to selection rules. 
The due population of the trion state during the pulse action 
is thus decreased. Therefore the contribution of the trion state 
to the expectation value of the Hamiltonian is lower than that 
for the no spin precession case. This leads to the discrepancy 
of dynamic phase. The reason for the discrepancy becom- 
ing more evident when Q > A is that more real population 
is transferred to the trion state during the pulse action. Thus 
the population transferred out of the two-level system results 
in more loss of contribution of the trion state to the dynamic 
phase. 
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FIG. 4. (Color online) Log-linear plot of dynamic phase component 
as a function of Q/A by two methods. The dashed and solid curves 
correspond to the result by method I and method II, respectively. 



IV. ELIMINATION OF DYNAMIC PHASES BY 
MULTIPULSE CONTROL 

Now we demonstrate the scheme to eliminate the dynamic 
phase by several In sec pulses control under slow precession 
approximation. Besides the symmetry of the dynamic phase 



to Q/A - +1, the overall phase and phase components are 
all odd functions of Q/A, as shown in Fig.|5] Utilizing these 
symmetry relations, we can choose particular parameters for 
each sec pulse, thus the dynamic phases induces by different 
pulses may cancel out. For the case of two pulses, we choose 
ri = Q/A for pulse 1, e.g. the red dot in Fig. |5] For pulse 
2, if r2 = -ri, the dynamic phase is the opposite number of 
that induced by pulse 1, a{-ri) = -a(ri) (e.g. the orange 
dot in Fig.|5]l. However, due to the odevity, the overall phase 
induced by pulse 2 is also the opposite number of that by pulse 
l,i.e. (fi{-r\) - -4>{r\),y{—r\) = —j{r\). Therefore after these 
two pulses all the phases cancel out. For the purpose of only 
eliminating the dynamic phases, we choose r2 - -l/^i = 
-A/Q for pulse 2. For this case the dynamic phase is still 
the opposite number of that induced by pulse 1 due to the 
symmetry (e.g. the green dot in Fig.|5]l, but the overall phases 
do not cancel out, leading to the geometric phase being added. 
Therefore the total geometric phase after two pulses reads 

Q A 
Jtoi = 01 -H 02 = 2 arctan(— ) + 2 arctan(-— ). (9) 
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FIG. 5. (Color online) Log-linear plot of dynamic and geometric 
phase components as a function of Q./A. The dashed (blue), dotted 
(purple) and solid (brown) curves correspond to the dynamic com- 
ponent a, overall phase (p and the geometric component y, respec- 
tively. Negative values of f2/A lead to negative overall and dynamic 
phases. The red, orange and green dots mark the values of dynamic 
phase when Q./A = 10, -10 and -0.1, respectively. Phases around 
Q./A = approach zero and are thus not sketched. 

Fig. |6] shows the total geometric phase y^, obtained after 
two pulses. It indicates that by our scheme the dynamic phase 
can be eliminated and arbitrary pure geometric phase can be 
obtained. Thus the rotation angle around z axis is a pure ge- 
ometric phase. In experiments, one may fix the pulse band- 
width Q, and only adjust the detuning from the trion reso- 
nance A to satisfy the parameter relation r2 = - 1 /ri for the 
two pulses. 

V. FIDELITY 

Fidelity is a direct measure to characterize how accurate 
a gate operation is implemented. The dominate mechanisms 
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FIG. 6. (Color online) The total geometric phase y^, obtained after 
two pulses, as a function of SI/ A. r = O/A, -A/Sl for pulse 1 and 2 
respectively. 

to deteriorate the fidelity of our nonadiabatic geometric rota- 
tion are the spin precession induced by the magnetic field and 
the population loss induced by the decay of the trion state. 
The former takes effect during the pulse action and between 
the two pulses while the latter functions only during the pulse 
action. The times of spin dephasing due to phonon and hy- 
perfine coupling are on the order of msec Il33ll and yusec fS?] 
respectively and thus could be neglected. Here we consider 
the effects of spin precession and population loss, and calcu- 
late the average fidelity by numerical solution of Schrodinger 
equation. The fidelity is defined as T(U) = K^|f/' f/,y|>F>p, 
where is the initial state, Uid and U are the ideal and actual 
operations respectively, and the average is taken over all input 
states [35]. We define / = U^Uid, same as Ref. lEIHl, then 
the fidelity is derived as 

where /,, , 7^, and Ijj are matrix elements of the operator I. 

In calculation, the actual operation U is firstly taken over 
the 3-dimensional Hilbert space then truncated in the 
subspace. Due to the effect of population loss, U is nonuni- 
tary. We take the parameters realistic in experiments as fol- 
lows: the magnetic field B - 0.29T [17] and the electron 
g factor \ge\ - 0.57 (4], which lead to the Larmar period 
T « 0.43ns, the pulse duration - 1.5ps ITtIi . and the trion 
lifetime t, - 900ps |37]. The two pulses are successively ap- 
plied and the distance between pulse centers is taken as 14t,/. 
We notice in Ref. llssll the square-wave pulse is adopted to 
realize geometric phase gates for QD charge qubits and the 
imperfection of the pulse shape is taken into account because 
it results in inaccurate pulse area. Here we consider the two 
pulses with a good hyperbolic secant shape and the pulse area 
is independent of the detuning ifTTll . In experiments, for the 
geometric rotation of n, only one pulse is required because it 
leads to a pure geometric phase, as demonstrated in Sec. III. 
The fidelities of other selected rotations of spin by multipulse 
control for studied QD parameters when B - 0.29T are listed 
in Table I. 



TABLE I. Fidelity of selected nonadiabatic geometric rotations about 
z for considered QD parameters when B = 0.29T. 



7wi (rad) 


Fidelity 


±7T/4 


99.38% 


±1x12 


99.43% 


±3;r/4 


99.48% 



The numerical results show the fidelities are equal for the 
positive and negative rotation angles which have the same ab- 
solute value. This is in full agreement with the analysis result 
in Eq. (|9]l: for +jtot only the sign of the detuning A is changed, 
the total stimulating effect of the two pulses are thus the same. 
The loss of fidelity is due to the combined action of the spin 
precession and population loss. Small angle corresponds to 
the detuning A close to Q for both pulses. For this case the ef- 
fect of spin precession is predominante due to the majority of 
virtual trion excitation. Large angle corresponds to one pulse 
with small A while the other pulse with large A. For this case 
the effect of population loss is enhanced (see Fig. |7]l due to 
the increased real trion excitation. Lower fidelity for smaller 
angle indicates that the effect of spin precession is stronger 
than population loss, which is more clearly shown in Fig. [8] 
When larger magnetic field is applied, the fidelity of small an- 
gle decreases more significantly than that of large angle. This 
indicates that compared with decay of trion state, the effect of 
spin precession is more detrimental to the fidelity, and weak 
magnetic field is a crucial condition for high fidelity of geo- 
metric rotation. 
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FIG. 7. (Color online) Total population after two pulses, as a function 
of the rotation angle for studied QDs. The transverse magnetic field 
is taken as S = 8T. 



VI. SUMMARY AND PROSPECTS 

In summary, we briefly review the fast spin rotation pro- 
posal proposed by Economou et al. [ 10] and analyze the ge- 
ometric and dynamic phase components of the overall phase 
induced by 2n sec pulses. The dependence of two phase com- 
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FIG. 8. (Color online) Fidelity of the operation as a function of the 
rotation angle for studied QDs. The upper (blue), medium (green) 
and lower (red) dots correspond to the fidelities when the transverse 
magnetic field is taken as S = 0.27T, 1.35T and 2.7T, respectively. 



without slow precession approximation. We discuss its phys- 
ical origin: the effect of spin precession. The spin precession 
leads to population loss out of the two-level system during the 
pulse action, resulting in the lower contribution of the trion 
state to the system energy. Utilizing the symmetry relations 
of the phases, we propose a scheme to eliminate the dynamic 
phase by multipulse control. By choosing the ratio Q/A as 
the negative reciprocal to each other for each pulse, the dy- 
namic phases induced by two pulse cancel out. Considering 
spin precession and population loss, high fidelity of the ge- 
ometric rotation is obtained when the weak magnetic field is 
applied. The detrimental effects of magnetic field and decay 
of trion state are compared and discussed. Our scheme real- 
izes high fidelity of geometric rotation of electron spin about 
an axis. Because of the evident advantage on short manipula- 
tion time and considerable feasibility of the proposal T IitIi . 
our scheme may generalize the proposal to universal nonadia- 
batic GQC. Further research may be focused on nonadiabatic 
two-qubit geometric gates in the quantum dot system. 



ponents on the ratio of the Rabi frequency to the detuning 
from the trion resonance is investigated. Numerical results 
indicate that if one pulse is applied, the overall phase is purely 
the geometric phase only when the spin is rotated by n. With 
other rotation angles the dynamic component is nonzero. An 
interesting effect is shown that the dynamic phase has a sym- 
metry as a function of rotation angle and Q/A. Results also 
show the discrepancy of dynamic phase for the cases with and 
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